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Abstract: We study the initial value problem for the wave equation and the ul- 
trahyperbolic equation for data posed on an initial hypersurface surface of mixed 
space- and timelike signature. We show that under a nonlocal constraint, the ini- 
tial value problem posed on codimension-one hypersurfaces — the Cauchy problem 
- has global unique solutions in the Sobolev spaces H m . Thus it is well-posed. 
In contrast, we show that the initial value problem on higher codimension hyper- 
surfaces is ill-posed due to failure of uniqueness, at least when specifying a finite 
number of derivatives of the data. This failure is in contrast to a uniqueness result 
for data given in an arbitrary neighborhood of such initial hypersurfaces, which 
Courant deduces from Asgeirsson's mean value theorem. We give a generalization 
of Courant's theorem which extends to a broader class of equations. The proofs 
use Fourier synthesis and the Holmgren- John uniqueness theorem. 
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1. Introduction 

The field equation 

Am - d 2 u = 

in Minkowski spacetime is of central physical importance, as it describes the prop- 
agation of many of the physical quantities described by field theories, including the 
components of the electromagnetic field in a vacuum. Its generalization to a theory 
which has multiple times is an ultrahyperbolic equation. The study of these equa- 
tions provides a useful window onto the mathematical status of physical theories 
involving multiple times, and perhaps more importantly, provides insight into the 
extent to which the ordinary concepts of causality and determinism survive the 
transition to multiple time dimensions. 

Consideration of theories with multiple times has been relatively rare because it 
is widely believed that they are inherently unstable, and thus are not deterministic 
in a physically meaningful sense. Certain significant developments in theoretical 
physics, notably string theory, require additional dimensions, and in most work to 
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datef the signature for the extra dimensions is spatial, reflecting in part this concern 
with instability. Motivated by this, the purpose of the present paper is to reconsider 
the questions of stability, uniqueness and determinism of the initial value problem 
in the presence of multiple time dimensions. We take the model field equation in 
this setting to be the simple generalizations of the wave equation to multiple times, 
the ultrahyperbolic equation. We find that the issue of stability and uniqueness for 
the Cauchy problem can be addressed by imposing nonlocal constraints that arise 
naturally from the field equations. 

It may be thought reasonable to go beyond the traditional Cauchy problem, and 
give initial data on hypersurfaces of higher codimension. We show that under the 
above constraints one can preserve stability in this setting, but uniqueness is lost, 
and thus determinism. Indeed, one may specify an arbitrary finite number of normal 
derivatives of the solution on the higher codimension hypersurface, and insist upon 
smooth solutions, yet still fail to achieve uniqueness. In contrast to this, we conclude 
with a result that essentially recovers and generalizes a theorem of Courant, which 
shows that the values of a solution in an arbitrarily small neighborhood of the initial 
hypersurface are sufficient to determine the solution uniquely In related and prior 
work, Woodhousc (1992) studied the case of two space and two time dimensions 
with initial data on a spacelike hypersurface (thus of codimension 2), using the 
Penrose twistor transform in the real case. He also recovered the uniqueness result 
of Courant and its implicit constraints on well-posed inital data for the Cauchy 
problem. Our work provides a rigorous analytic alternative for his solution method, 
which is not restricted to this choice of space and time dimensions. We remark that 
none of these results rely upon properties of analyticity of the data or the solution. 

To fix our notation, the wave equation in di-many space dimensions and one 
time dimension is 

di 

A x u-d 2 y u:=J2d 2 x .u-d 2 y u = . (1.1) 
j=i 

The standard Cauchy problem is posed on N = {(x, y) G M^ 1 x WL y : y = 0}, a 
spacelikc codimension one linear hypersurface, for initial data 

u(x,0) = f(x), d y u(x,0) = g(x). 

A nonstandard Cauchy problem is posed for a linear hypersurface of mixed signature 
N = {(x, y) : xi = 0} C x R y , namely 

w(0, x', y) = f(x', y), d Xl u(0, x', y) = g(x', y), 

where the notation is that x = (x\,x') € M. dl . Courant (1962) calls this the 
non-spacelike Cauchy problem, but to avoid confusion with the non-characteristic 
Cauchy problem, we call it a Cauchy problem of mixed signature. 
An ultrahyperbolic equation has the form 

d% d 2 

A x u-A y u:=J^d 2 .u-J^d 2 y .u = 0, (1.2) 

f Exceptions include the work of Tegmark (1997), Hull (1999), Hull & Khuri (2000), and Bars 
(2001). 



On determinism and well-posedness in multiple time dimensions 



3 



where x € R dl are considered to be the spacelike variables and y <G R d2 are time- 
like. The Cauchy problem considers initial data posed on a linear hypersurface of 
codimcnsion one. Choosing yi as the direction of evolution, Cauchy data consist of 



on the hypersurface N = {(x,y) € R dl x R^ 2 : y x = 0}. 

The initial value problem on a higher codimension hypersurface M could take 
various forms. A natural problem from the perspective of theories with multiple 
times is to consider the spacelike hypersurface M = {(x, y) £ Rf. 1 x R d2 : y = 0} 
of codimcnsion d 2 . Alternatively, one may consider more general M = {(x,y) € 
M^ 1 x R d2 : x Pl+1 = ■■■ = x dl = 0, yp 2+ i = ■■■ = y d2 _ 1 = 0} where < pi < d x 
and < P2 < d 2 — 1. There is in either case a question as to how much data, and 
what constraints, are to be considered on M. Some of the options are to (i) give 
the zeroth and first normal derivatives of w on M, (ii) give some finite number of 
derivatives of u on M which are compatible with the constraint imposed by the 
ultrahypcrbolic equation, or (iii) specify infinitely many derivatives of u on M. In 
this paper we consider the first two of these cases. 

An outline of the results of this paper is as follows. In section 2 we use Fourier 
methods to show that the Cauchy problem for the ultrahyperbolic equation (1.2) 
is ill-posed in general but well-posed on Sobolev spaces H m if an explicit nonlocal 
constraint is imposed upon the Cauchy data. This applies as well to the wave equa- 
tion with Cauchy data on a mixed signature hypersurface. In section 3 we consider 
the initial value problem for data given on higher codimension hypersurfaces, and 
we find that solutions are highly nonunique for the initial value problems of type 
(i) and (ii) above, even among H m smooth solutions and with the imposition of the 
constraint given in section 2. In particular, for theories with multiple times that 
can be transformed to the form of equation (1.2), data posed on the hypersurface 
M = {y = 0} do not uniquely determine the solution at any other point in time 
y £ R d2 \{0}. The extension problem for higher numbers of derivatives is treated by 
the same method as case (i) of zeroth and first normal derivatives. Regarding case 
(iii), in which one specifies infinitely many derivatives on the initial hypersurface 
M, we do not have an answer. We do show in section 4 that among smooth solu- 
tions, data in an arbitrarily small ellipsoidal neighborhood of a disk in M uniquely 
determine the data in the envelope of its light cones. This is analogous to a result 
in Courant (1962) that is derived from Asgeirsson's mean value theorem. 



Let x G R dl and y e R d2 be the Cartesian coordinates of space-time, denote 
y = (yi,y') and consider the Cauchy problem of evolution in the coordinate y\. 
The Cauchy problem of mixed signature that we address is posed as 



with Cauchy data u{x, 0, y') = u (x, y') and d Vl u(x, 0, y') = ui(x, y'). The standard 
Sobolev spaces H m of functions of the variables (x,y') are defined as closures of 



u(x, 0, y') = f(x, y'), d Vl u(x, 0, y') = g(x, y') 



2. The Cauchy problem 




(2.1) 
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C§°(R dl x M d2_1 ) with respect to the norms 

ll/t= E f\d^,f(x,y')\ 2 dxdy' . 

\a\ + \/3\<m J 

Additionally, there is an energy functional, of indefinite sign, that is associated with 
equation (2.1), namely 



E{u) := -J J \d yi u\ + \V x u\ - \V y 'u\ dxdy' . 

Theorem 2.1. Suppose that the evolution mapping y\ — > (g"^' 1 ^^) * s inC 1 ( 
H 1 x TTien £/ie energy is conserved along a solution u(-,yi, ■): 

£(«(■, j/i, ■))=£(«(■,(),■)). 

Proof. Given L*' 1 ; 91 ' 9 e C 1 , the following calculation is justified: 

\ Oy ^ u y x , y i , y )/ 

d yi E(u) = J[(dy 1 u-d% 1 u + V x u-V x d Vl u-'Vy>u-'Vy>dy 1 u)dxdy' 
= II dy^d^u + A x u + A y ,u)dxdy' 



□ 



The key issue is that the Cauchy problem above for equation (2.1) is ill-posed 
for di > 2 and solutions are not in general in C 1 (M. yi : H 1 x H°). The energy 
is indefinite and in particular not bounded below, hence it does not in general 
define an energy norm with which to control the Sobolev norms of solutions of the 
evolution equations. 

To move to the next level of analysis, we give a Fourier synthesis of the evolution 
operator for the Cauchy problem of mixed signature. Given ("°) e x H m , 

consider the Fourier space variables (x,y') — > (£,7]') and define the Fourier trans- 
form in the standard way, 

m^v')\_ i rr e -i^ e -i,'y'( u oM)\ dxdy/ 



where d = {d\ + di — 1). On a formal level equation (2.1) under Fourier transform 
will read 

d 2 yi u=(-\Z\ 2 + \r 1 f)u, 



giving rise to the expression for the propagator, exp(yi^/ A x — A y >). The solution 
thus reads 

, u , / lf|2 , ,,2 ^ Ic , M\f\t\ 2 ~W\ 2 yi) ^ 

u{£,yi,V ) = cos(y |£| -\t)'\ yi)uo{£,r)) + uxi^rj) 

^\tf-w\ 2 
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for \r]'\ < while 



; T sinh(V Ir/I 2 — \£ \ 2 y\) 

vVl 2 -l£l 2 



for |£| < \rj'\. That is, the dispersion relation 



u&rf) = \/lC| 2 - W\ 2 (2-2) 

holds in the Fourier space region {\rf\ < |£|}, while in the complementary region 
the evolution of a Fourier mode is described by the Lyapunov exponent 



mv') = vVi 2 - iei 2 . (2.3) 

When the propagator is applied to data ("°) which is analytic, this solution 
exists for at least short time; for analytic data of exponential type, the solution is 
global. However, it is clear that general initial data in H m+1 x H m do not even 
give rise to solutions which are tempered distributions for any nonzero y\. 

On the other hand, imposing a constraint on the initial data, the solution process 
is well defined. The fact that some constraint is necessary is indeed evident from 
the Asgeirsson mean value theorem, and its consequences, as discussed in Courant 
(1962). The form of this nonlocal constraint is evident from the Fourier synthesis, 
as we shall now see. 

Define a phase space X using an energy norm adapted to the propagator of 
equation (2.1). Using the definition of the dispersion relation (2.2) and the Lyapunov 
exponent (2.3), and the Plancherel identity, set v — and 



\v" 2 



:= J J W 2 (£,r/)|t>o(£,r/)| 2 <« 



x 

! </'l ; ! 

A 2 (£,r/)M£,?/)| 2 « 

{\t\<\v'\} 



This is a norm, unlike the actual energy associated with the equation (2.1), and can 
be used to control solutions when the propagator is restricted to the appropriate 
stable and/or unstable subspaccs of X. Define 

* s ={*= (:) - x -\ + Wy) = ° for |c| < wi ) (2 - 4) 

X u ={v€X: I^ o( ^r ? ')-||^)=Ofor |£| < |r/|} (2.5) 
X C - {«6 X : supp(!°) %rf) C > |r/|}} 

= x s nx u . 



and 
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The subspace X s corresponds to the center stable subspace for evolution in y\ £ 
R+, the subspace X u corresponds to the center unstable subspace, and X c is the 
center subspace. This nomenclature is supported by the following theorem. 

Theorem 2.2. For ("°) £ X s , the Cauchy problem of mixed signature for equation 
(2.1) has a unique solution in X for all y\ £ M + . For (^°) £ X u the problem has 
a unique solution for all y\ £ R~ , and whenever (^°) £ X c the solution exists 
globally in y\ £ K. In each of these cases, the map y\ — > u(x,y\,y') is C 1 . 

Denote the propagators by <1> S , $ u and & c for data in the respective subspaces. 
These solutions are continuous with respect to their Cauchy data taken in the 
respective subspaces. This is the result of the next theorem. 

Theorem 2.3. Given two phase space points u = ("°),t> = ("°) £ X s , then for 

yi >o, 

II I 

The analogous estimate holds for u,v £ X u , for yi < 0: 



I 2 II II 2 



(2.6) 



K>)-<>)|| 2 <||«-«" 2 



X 



(2.7) 



For u £ X c , $^ = $^ for y x > and $^ = <j>^ /or j/i < 7 and equality holds in 



^ Vl 



Vl 



-yi 



both (2.6) and (2.7). 

Proof. It suffices in Theorem 2.3 to prove the first statement. In X s the solution 
has two components, distinguished by their Fourier support. Consider first ("°) 
such that supp(uo,Mi) C {|£| > \n'\} := R\, which gives the center component of 
the evolution. The propagator is expressed 

sin(t^yi) 



Vl 



COs(0J7/i) 

— oj sin(wj/i) cos(wyi) 



where w 
norm, 



w(£, r/) is the dispersion relation (2.2). Evaluating this in the energy 



sin(wyi) „ 
cos(w2/i)uo H ui 



+ |— iv sin(ujyi)uo + cos(uiyi)ui\ d^dr]' 

= JJ{\u \ 2 U 2 + \u 1 \ 2 )d^d7 1 ' 



(2.8) 



x 



The propagator on the complementary space is more sensitive. Let us suppose that 
supp(Mo,Mi) C {\n'\ > |£|}, then A(£,r/) > and we express the propagator in 
terms of its Fourier transform as 



cosh(A yi ) sinh [ Xyi) 
Asinh(Ayi) cosh(Ayi) 



A 1 



5 -Ai/i 



1 -* 

-A 1 



Uq 

Ml . 



On determinism and well-posedness in multiple time dimensions 



7 



The subspace X s consists of precisely those data which lie in the null space of the 
first term; this is the expression of the constraint 



A«o(£,f/) + «i (£,»/) = 
Measuring the remaining term in energy norm, we find 

-2\ yi 



(2.9) 



s ( u o 

Ul 



X 



uo-- 



A + l-Auo + mi| 



//■ 

J J e- 2 ^ (|u | 2 A 2 + |tii| 2 ) dZdrf . 



d^dn' 



Since we consider the propagator $^ for j/i > 0, the exponent — \y\ is negative, 
and therefore 



<f>' 



< 



x 



for tt 



€ X s . For general data in X s , one decomposes it into its components 



with support in {|£| > |r/|} for which we use (2.8), and its component supported in 
{Wl > |£|}> which in addition satisfies the constraint (2.9). Therefore on X s 



|*£(«)|£<N& 



Bounded operators on X s arc continuous with respect to u e X 5 , and it is easy 
to see that the solution behaves continuously with respect to y\ > as well. The 
case for the subspace X u is proved by the same arguments, after reversing time 
yi — > —yi. This proves Theorems 2.2 and 2.3. We remark that on the center 
subspace X c , which yields global solutions, both constraints are imposed 



implying that u (£,f?') = = fii (£,»/) for all < |tj'|}. 



(2.10) 
□ 



The proof extends to the initial value problem posed in higher energy spaces, 
defined by 



IMI 



X™ ■' 



\a\ + \P\<m 



X 



We then have 



Corollary 2.4. The higher energy space X m decomposes into three subspaces, 
X m - S , X m - U and X m ' c = X m ' s n X m - U such that for it, v € X m ' S and y x > 



\$ s (u) - $ S (v)\\ < 
I yi \ U J yi^ u > \ \ x m — 



u — v\ 



X'' 



while for y\ < and u, v e X m,,7 ; 



< lit 



For u,v e X m - C both estimates hold, and a global solution exists which has prop- 
erties of higher Sobolev regularity. When m > ((d\ + (d 2 — l))/2) + 2 then such 
solutions are known to be classical C 2 solutions by the Sobolev embedding theorem. 
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It is natural to estimate solutions with respect to the energy norm; indeed, it is 
the energy when restricted to the center subspace X c . Thus the problem is well- 
posed in the following sense: data in X s continuously propagates to all yi € R + , 
data in X u continuously propagates to all times yi e R~, and data in X c \ which 
constitute an infinite-dimensional Hilbert space, are defined globally in time. In the 
case of the ordinary wave equation (d\ = 1), solutions in X c correspond to the full 
energy space H 1 x L 2 . 

3. The initial value problem in higher codimension 

In the presence of multiple time dimensions, spacelike hypersurfaces are necessarily 
of higher codimension. Therefore one might consider the initial value problem with 
data posed on a hypersurface of codimension greater than or equal to two. Such 
problems are generally ill-posed. Indeed, we show that solutions can be singular 
for standard classes of data. Moreover, even imposing the constraint discussed in 
section 2, which is the requirement of global existence, smooth solutions are highly 
non-unique. The purpose of this section is to study the extension problem of data 
posed on a non-degenerate higher codimension hypersurface M to Cauchy data 
on a codimension one hypersurface N. There is a lot of freedom in choosing this 
extension, even under the constraint equation (2.9) on the resulting Cauchy data. 
Other extensions can be chosen to fail to satisfy the constraint. Thus the initial 
value problem fails to be well-posed in several ways: resulting solutions may be 
singular, or they may be selected to satisfy the constraint and be regular for all 
yi e R, however they will not be unique. 



Our analysis is illustrated in the example case of M = {y\ = y 2 = 0} and 
N = {yi = 0} subspaces of R 3 . We suppose that initial data for a solution u(x, y) 
is given on M in the form 



where wq{x\) — u(xi,0), wi (xi) — d Vl u(xi,0) and wqi(xi) — d V2 u(xi,0), corre- 
sponding to the values of the solution and its normal derivatives on M. The object 
is to extend w(x\) to Cauchy data (uq{xi, y%), Ui(x\, 1/2) on N which satisfies 



We give extensions which satisfy the constraint (2.9), therefore giving rise to global 
solutions in y\ £ R. Such extensions are nonuniquc. Additionally, there are exten- 
sions which fail to satisfy the constraint, lying in X S \X C or X U \X C or neither. 

Definition 3.1. The extension operator is given by 



(a) Codimension 2 to codimension 1 in R 3 



w{x\) = (w (x 1 ),w 10 (x 1 ),w 01 (x 1 )) 



u (a;i,0) = w {xi) 
ui(xi,0) = tuio(a;i) 



and the compatibility condition 



<9j, 2 «o(zi,0) = woi(zi)- 
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where the kernel function rj') is chosen such that for all 

1 S xkU)drj 



In order to satisfy the constraint, we ask additionally that supp(x(£, ?/)) 5= {W\ < 
|£| }. A reasonable choice is to take 

xM) = ^'/l£l)^, 



/or V(6>) € C^°([-l, 1]), V(6>) > even, and 



(3.1) 



Theorem 3.2. The extension operator E is a bounded operator on the following 
space of functions: 

E : £T _1 / 2 (M) -» L 2 (7V) 

(iT 1 / 2 n H m - l ' 2 ){M) -» iT"(A0 . 

/n addition, when w e H~ 3 / 2 (M) then y2E(w) e L 2 (N) and furthermore 

y 2 E : H m ~ s ' 2 {M) -» iT™(iV). 

Using the extension operator, we generate constraint-satisfying Cauchy data on 
iV from initial data on M as follows: 

uo(a;i,2/2) := E(wa){xi,y 2 ) + y2E(w i)(xi, t/ 2 ) 
ui(a;i,2/2) := ^(wio)^!,^)- 

Checking that this is a legitimate choice, we have 

«i(xi,0) = ^// e* B1 «>io(Ox(£,'/)<W 

= wio(xi) 

because of the normalization (eqn 3.1) of ^. Similarly, 

u o (a;i,0) = E(w o ){xi,0) = w {xi). 

The compatibility condition is satisfied, since 

d y2 uo(xi,0) = d y2 E(w Q )(x 1 ,0) + E(w )(x 1 , 0) 
= d V2 E(w )(x-i,0) +w i(x 1 ). 
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The first term of the RHS vanishes because 



1 



/2tt 
0, 



where we have used that J 9ip(9)d9 — because ip(-) has been chosen to be even. 

The pair of functions (uo(xi, 1/2), ui(xi, 1/2)) gives Cauchy data for the codi- 
mension one problem that is discussed in Section 2. Because of the properties of 
the extension, it satisfies the constraint conditions of X c for solutions which are 
globally defined in y\. In order to apply the existence theorem, the energy norm of 
this Cauchy data must be finite. 

Theorem 3.3. Suppose that w G H X I 2 {M), w 01 € H~ 1/2 and w 10 G H^ 1 / 2 . 
Then the energy norm of the extension uq = E(wo)(xi,y2) + ?/2-E(?«oi)(:ei, 2/2), 
ui{x 1 ,y 2 ) = E(w 10 )(x 1 ,y 2 ) is finite: 

\\(u ,Ui)\\ 2 x c < C(\\w \\% 1/2 + \\w l\\ 2 H-i/2 + HwioH^-1/2). 

Additionally, the higher energy norms with which one defines the X m topology for 
(u ,ui) are also bounded by this extension process, namely 

\\(U0,Ul)\\ 2 xm ,c < C m (\\w \\ 2 Hm+l/2 + ||W0l||^ m -l/2 + ||wio||^™-l/2). 

Proof, (of Theorem 3.2): Using the Plancherel identity, the L 2 (N) norm of E(w) is 
\\E{w)\\ 2 L , (N) = J J \u,{0\ 2 ^{v'/\t\)^dn'dt 

= J ^N£)I 2 (/^V/I£I)^')^ 

since 9 = rj j |£| and 

1 

/ V>V/l£l)^/ = J 4, 2 (0)de. 

-1 

The identity extends to the Sobolev space H m (N) ; it suffices to calculate 1 1 9™ E(w) 1 1 L2 
and \\d™E(w)\\ L2 : 

\\d™E( W )\\ L2{N) = J J K ? )| 2 |e| 2m V 2 (Wiei)^2« 

= / HOI^r- 1 (/ V>V/KI)^*/) di 

L 2 [-l,l] \\ w \\h^- 1 / 2 (M) ■ 
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The second quantity is similar: 

\K E H\\ LHN) = JJ KOI 2 KI 2ro V<V/ KDji W\ 2m dr,'dt 



2m 



= C m \\w\ 



ffm-l/2( M ) , 



where C m = / 2m ip 2 {0)d9. The third and fourth estimates of the theorem involve 

-l 

y 2 E(w), whose Fourier transform is 
Measuring the L 2 norm of y 2 E(w), 



\\V2E(w) 



L 2 (N) 



fB^'l\li\)— 2 



dn'di 



// 



H- :, / 2 (M) 



with C = / | t/^ | 2 dO. The calculations of the H m norms of y2E(w) are similar. □ 
-l 

Proof, (of Theorem 3.3): Given initial data (wo,Woi, Wio){xi) we are to give con- 
ditions under which the energy norm of the extension (uo,ui) is finite. First of all, 

111 1 1 2 

the contribution to the energy given by u\ is simply ± ||ui|| i 2, hence by Theorem 
3.2 it is bounded by C ||i«io||jj-i/2- There are two contributions from uq, which can 
be expressed using the Plancherel identity: 

JJu; 2 (t,V , )\MO\ 2 X 2 (Lv , )dv'dZ + JJw\^ „') |u, 01 (£)| 2 |^ X 2 (£, //) | 2 « • 
Using that x(€>v') — i'iv'/ > we estimate these two integrals: 

/ f (iCI 2 - Wf) MOI Vfa'/ IZD-^dn'dt 



W\<\t\} 



\M0\' 



<C\\w \\% 1/2 . 



ki-^Wv/ki)^/ 



dt 
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J f (iCI 2 - \vf) Ni(OI 2 tfV/ 



{|VKI«|} 




< C|koi||jy-i/2 • 



□ 



(b) The extension problem for general spacelike data 

We now consider the general problem of initial data given on a maximal spacelikc 
hypersurface of dimension d\, extending it to Cauchy data on a codimension one 
hypersurface. That is, for (x,y) £ R^ 1 x 



Initial data on M take the form w(x) — (wo(x),w a (x)) where a solution u(x,y) of 
the field equation (1.2) is asked to satisfy 

u(x,y) = w (x) 

with its first derivatives normal to M satisfying 

dyu(x, 0) = w a (x) 

where a e N d2 is the multi-index a — (a\, ad 2 ), \a\ — 1, such that only one 
a.j = 1 and the rest are zero. The object is to extend w(x) to Cauchy data on N 
while satisfying the constraints (2.9) to be in X° . This Cauchy data satisfies 

uo(a;,0) = w (x) 
u a >(x,0) = W0a'(x) 

for a' = («2, ctd 2 ) and the first derivatives normal to N satisfy 

d Vl u(x,0) = wi (x) . 

Following the construction given in section 3.1, define an extension operator 



To satisfy the constraint that E(w) G X c , we ask that supp(x(^, n')) C {(£,r)') : 
\v'\ < l£|}- Such kernel functions are readily constructed (they are far from being 
uniquely determined). For example, a variant of our construction of section 3.1 is 



M = {y = 0} C N = {y x = 0}. 




where the kernel function is even in 77 and satisfies 
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based on choice of a Cq 3 function ip(9) > 0, with supp(^) C £>i(0), the ball of 
radius one. Then define 

x(Z,v') = i>(v'/\t;\)j^. 

We note that \ lSi even m V if V'W is even, and that 




This is normalized to one by choice of rp. 

Theorem 3.4. The extension operator E is bounded on the following function 
spaces: 

E : H^(M) -> L 2 (N) 

H^(M) n H m+l ^{M) -» ff m (iV) , 

m the exponent of Sobolev regularity, and 

y'E : H^^iM) -» L 2 (A0 

H z ^ A (M)nH m - 1± ^{M) -> H m {N) . 

Using the extension operator the vector function w{x) — (w (x) 7 w a (x)) 
extends to Cauchy data on N as follows: 

u Q {x,y') := E(w )(x,y') + ^ y a ' ■ E(w 0a ,){x, y') 

\a'\=l 

ui(x,y') := E(w 10 )(x,y') . 

This is seen to extend the initial data w(x) in the required way, and in addition 
it satisfies the constraint that (uo,u\) G X c . However, measuring the functions 
(uo, tti) in the energy norm is more appropriate for the Cauchy problem, hence we 
also state estimates in this setting. 

Theorem 3.5. Given w G H(M) and w a G H(M), the energy norm of the 
extension 

(u , ui) = {E(w ) + y a ' ■ E(w(p, a ')),E(w w )) 

is finite; indeed 

\\(u ,U 1 )\\ 2 x c < C(||uio|| 2 3-d2 + ||w0a'|| 2 i=*2 + ll w lo|| 2 i-Ja ) • 
H 2 H 2 H 2 

The proofs of Theorems 3.4 and 3.5 are similar to the proofs of Theorems 3.2 
and 3.3, to which we refer the reader. 
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(c) The extension problem for mixed spacelike and timelike data 

As a final case, we consider the extension problem for initial data on a lower 
dimensional hypersurface M of mixed signature. Given zero'th and first normal 
derivatives of a solution u(x, y) on M, the object is to extend this data to the 
codimcnsion one hypersurface N — {y\ = 0} in such a way that the constraint 
for well-posedness is satisfied. This is not always possible for arbitrary data w — 
(wo,w a ) posed on M, due to analogous lower dimensional constraints on M. But 
it is possible, along with attendant Sobolev bounds on the extended functions, in 
most cases. This situation will be explained below. 

To set the notation, we consider spacelike and timelike coordinates on M to be 
(x , y) e MP 1 x W 2 , with their Fourier transform variables denoted (|, ff) e W 1 x W 2 . 
The complementary variables will be denoted (x",y") G x R d 2-P2-i an( j 

(£",V") G K dl_Pl x R*-P2-i 5 so that coordinates on N arc (x,y') = (x,x",y,y"). 
The evolution variable remains y\ . 

Initial data for a solution u(x,y) is given on N, which is expressed in the form 

(u,d^'/,u,d^u,d^,u)(x,y,0,0) = (w ,w a „,W(3 1 ,w f )„)(x,y),wherea" = (a Pl+1 , ...,a dl ), 
f3" = (/3 P2+ i, Jid 2 ) are multi-indices such that |a"| + |/3"| + |/3i| = 1. The idea is the 
same as in sections 3.1 and 3.2, namely to extend (wo, w a », , wp») to constraint- 
satisfying Cauchy data on N in such a way that a solution u(x, y) — u(x, x" , y, y") 
to the field equation (1.2) satisfies 

u{x,0,0,y, 0) = w {x,y) 

and 

d Vl u(x,0,0,y,0) = w o01 (x,y) , 
as well as the compatibility conditions 

d"",u(x,0,0,y,0) = w (a „ fi) (x,y) 
dy„ u(x, 0, 0, y, 0) = w ( o,/3") (x, y) 

The existence of such an extension follows as in Theorems 3.4 and 3.5 from the 
construction of an extension operator E with certain boundedness properties on 
appropriate Sobolev spaces. We will focus our analysis therefore on the extension 
operators. 

Again following section 3.1, define an extension operator 

E(w)(x,y') = / / x(lC,V,v")^"d V " = 1. 

Furthermore, to satisfy the constraint that E(w) € X c for arbitrary data w, we 
ask that 

supp( X (e,V))c{(e,r;'):|r/'| 2 <|e| 2 }:=i?i. 

These two conditions are always satisfiable, except in the case £" = {0}, meaning 
that d\ = pi and the extension subspace {(£,", rj")} is purely timelike. 
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It is to be expected that the constraint induces a restriction on the data w(£, 77) 
in the vicinity of the "lightcone" ||<f| = C M. Subdivide M into two sets, 

R, :={{lf,)&M:\f,\<\z\) 
R 2 :={(|^)eM:|^|>|||}. 

The orthogonal projections onto functions supported in R\, R2 respectively, are 
denoted 7Ti and 7r 2 . We use standard Sobolev spaces to quantify data supported in 
Ri, namely 



H r = { 



w(x,y) G rangc(-7ri) : \\w 



2 „ 2 



Mtv) U + \v\ ) n dtdTj < +co 



Hi 



Over i? 2 we use a modified form of Sobolev norm which is given by 



w(x,y) £ range(7r 2 ) : \\w\\ 2 K r = J J w{£,v) 



(I +m 2 r - 

2 d^dri < +00 



H 2 - 



where 

e := d 1 + d 2 - (pi + p 2 ) - 1- 
We note that in the case where di = pi, #1 = {0} and K n = £r-§(d 2 -i>2-i) . More 
generally, define 



K = { 



w(x,y) G rangc(7r 2 ) : \\w 



77) 



2 u 2 +w 2 r 



fl 2 



l^l 2 



-rf^dry < +00 



Decompose an arbitrary function w — ttiw + tt 2 w, so that its components possess 
Fourier support in Ri and i? 2 respectively. 

Theorem 3.6. If di > p\ then there is a choice of kernel x (indeed there are many 
such choices) such that u — E(w) satisfies 

\\uf L 2 < C(lkiw||^-^(e 0) + ||7T 2 w||^ ). 

Higher Sobolev norms of ' u — E(w) are bounded as follows 

\\uf H r < Cr(lkiu;||^ r _i (eo , + h2w\\ 2 Kr . 

In case d\ = pi, it is not possible to extend arbitrary data to a function u — E{w) 
which satisfies the constraint supp(u(£, 77')) C R 1 . However, if initially supp(u)(£, 77')) 
Ri (i.e., w = 7T2W), then such an extension is possible, and we have, for u = E{w), 

\\uf L2 < C\\wf K0 , 



|| U ||* < C r \\W" 2 



K v 
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Proof. The proof of Theorem 3.6 depends upon the construction of a kernel x(£, fj, 77") 
with satisfactory properties. This construction is slightly different in the two differ- 
ent regions of Fourier space 

Ri ■= {(Iv) ■ \v\ < |l|} and R 2 := [{If,) : \fj\ > |||} 

where we note that the region R 2 contains the data which would lead to an ill-posed 
initial value problem if M were considered itself as a codimension one hypersurface. 
To extend data posed on region Ri, define 



r," 



(I 



U + \v\ 2 ) 



7) 



1 



(I 



+ 1^ 



where tpi(0\, 62) is a Qf 3 function of (rfi — pi) x (d 2 — p 2 — 1) variables, respectively, 
with support in the set 1 6*2 1 < |#i|. Therefore Xi has support in the set 



implying that 



\V 



(i 2 + w 2 )? 



I -|2 



> 



~ 2 9 1 



< c 



_ 2 



+ m 2 = ier 



This is the appropriate region of support from functions v = E(w) to lie in the 
constraint-satisfying subspace of L 2 (N). In order that E be an extension operator, 
we furthermore require that 



2tt 



l+d 2 



'-II 



Xi(i,V,e,v")dt"dri" 



// 



l^l 2 )^ 



+ l^l 2 ) 



7) 



\fj\ 2 )^ ea 



= // M0i,0 2 )d9 1 d6 2 . 



l+d 2 



Asking that this latter integral equal the normalizing constant V27T , and 

asking for ipi to be even in its variables (61,62) gives an acceptable kernel for the 
extension operator. We note again that this choice of kernel is highly nonuniquc. 



}■ 



we can also attempt a construction 



On the region R 2 = {(£, fj) e M : \ fj\ > 

of our extension operator. By itself, this region would give rise to data in L?(M) for 
which the Cauchy problem of mixed type is ill-posed. The extension operator will 
nonetheless come up with data u — E(w) for which the well-posedness constraint is 
satisfied, if this is possible. That is, as long as di > pi, so that {£"} is not restricted 
to the zero-dimensional vector space, extensions can be found in a way that the 



< I fj\ can be made 



default in satisfying the constraint caused by the fact that 

up with a choice of large In practice, we will build X2(£, V, £"> v") so t na t its 
support is in the regions 

{|»7l>|l|}=-R2 
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as well as 



, i„//|2 



+ wr< c +ri ; 



implying that < (\fj\ 



) + (I 7 ?" 1 2 + l£"| 2 )- Thus we require c?i > pi. Following 



the above examples, assume that d\ > p\ and set 



for (£,77) <E i?2- Let ip 2 {0i,9 2 ) be a Cg function of e = rfi + d 2 — {p\ + P2) — 1 
variables, as before and require that 



M0u02)d6 1 d9 2 



r 2 , J 



1 



l^l 2 



Furthermore, ask that ip(9\, 2 ) be even in (61, 9 2 ). Finally ask that the support of 
^{0i,0 2 ) be in the set 

{(0i,0 2 ) -el -el > i}. 

Such requirements are satisfied by many possible choices of ij). In doing so, we 
arrive at a satisfactory kernel of an extension operator E with the property that 
all functions u = E(w) in its range have Fourier support satisfying supp(u) < 

1 1 T]' 1 2 < I £| 2 j. The singularities introduced at the boundaries of the lightconc j \rj\ = 

M by the kernel \i impose more severe constraints on the functions w that are per- 
mitted in the domain of the operator E; this is the origin of the somewhat unusual 
requirements on functions w(x,y) from which we can reasonably draw our data. 
The Sobolev estimates of the proof are similar to those of Theorems 3.4 and 3.5 
and we leave the details to the reader. □ 



Finally, we show that a sufficiently large class of data (w , W( a " fl), w {o,p")i w i3i) 
on M extends to Cauchy data on the hypersurface N which is both of finite energy 
and satisfies the constraint. This extension is given by 

u (x,y') = E(w )(x,y')+ ^ x" a " E(w {a „ fi) )(x, y') + ]T y" " E^p^x, y') 

l«"|=i I/?" l=i 

(3.2) 

ui{x,y') = E{w {0 ^ l) ){x,y l ). 

By design, this Cauchy data satisfies the constraint, that is, (uo, u i) <= X c , the cen- 
ter manifold. As before, its restriction to M reduces to the data (wo, W( Q ",o) > w {o,i3")){x 
The only remaining task is to show that its energy norm is finite. Recall that in 
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this context the energy norm is 

H{uq,Ux) = i j j \ui\ 2 + |V x u | 2 - \\7 y ,u \ 2 dxdy' 

N 

= \JJ \Mtv')\ 2 + m 2 -W\ 2 )Mtv')d(d v '. 

N 

To show that this energy is finite for the extension (3.2), we use the results of 
Theorem 3.6. 

Theorem 3.7. Given data (w , W( q ".o)j w (o,,3"): w /3i) on M with \a"\ — \(3"\ = 1, 
suppose that 

hiW \\ H e +i + ^2 || 7r i U7 (a",o)|| H e 0+ i + X! IK^o./nlli^o+i < + 00 ( 3 - 3 ) 

l«"|=l 1/3" 1 = 1 

h2Wo\\ K i + W^Wia" ,0)\\ K i + IK 2U; (0,/3")llxi < +°° ( 3 ' 4 ) 

l«"| = l |/3"|=1 

and 

Iki^ftll^o + ll^^ftll^o + oo. 

TTien i/ie extension (uq,Ui) given by expression (3.2) has finite energy and lies in 
the center subspace X . If di — p\, then we have to ask that iT2W 7 — in the above 
statement, for all multi-indices 7 in question. 

Proof. Estimates on the contributions of wq to uq follow immediately from Theorem 
3.6, as do the estimates for u\ = E(wp 1 ). Therefore we only have to consider 
contributions in one of the two possible forms: 

x" a "E(w {a ,, fi) ) \a"\ = 1 



or 



y"P" E{w {Ot0ll) ) 



1 . 



The energy norm includes the quantities 



and 



y"0"E{w {Ot0ll) ) 



since the estimates are similar we will give a sketch of one of them. 



H 1 



< 



x" a "E(w {a „ fi) ) 



2 

H 1 



-dfEiw^^mf + wiy. 



_ 2 



_ 2 



L 2 
1 



- 2 



+ 



f" 77" \ 1 12 



7riu;( a »,o)(|, 77) 



I ~|2 



7) 



I ~|2 



T2^(a'',0)(^.»?) 



2 _ 



d£drjd?'drf' 



The ^"-derivative introduces one extra factor of ( £ + 1 77 | 2 ) , respectively ( | 77 | 2 - 



), 



into the denominator. The integral over (£", 77") gives a constant, depending upon ipi 
and 7^2, as a bound, while the resulting integral over the variable (£, 77) is bounded by 
the iJ 1_e ° norm (respectively, the K\ norm) of ui( a n y This finishes the proof. □ 



On determinism and well-posedness in multiple time dimensions 



19 



4. Failure of uniqueness in higher codimension 

The question addressed in this section is the uniqueness of solutions with prescribed 
initial data on a hypersurface M of codimension greater than one. This is a nontriv- 
ial issue if one requires that solutions exist globally in space-time, which has been 
the focus of the analysis in the preceding sections. In section 3 we showed that initial 
data consisting of the values of the solution u(x, y) and its first normal derivatives 
on M, through a procedure of extension, give rise to constraint-satisfying Cauchy 
data on a codimcnsion-one hypersurface N . These extensions are highly nonunique, 
and therefore so are the resulting global solutions. 

We now raise the question whether prescribing an arbitrarily large but finite 
number of normal derivatives on M, as well as insisting upon global solutions, would 
remedy the nonuniqueness. This data should satisfy the compatibility conditions 
implied by the commuting of mixed partial derivatives and by equation (1.2). Given 
Courant's classic result (1962) in the case of purely timclike M, that data given in 
any e-tubular neighborhood of M within N determine solutions uniquely in the C 2 
category, one might think that specifying additional data for u(x, y) on M would 
suffice. In fact, if one specifies any finite number of derivatives of u on M it does 
not. 

Theorem 4.1. Given k, there exist constraint-satisfying data Uo,U\ on N which 
vanish to order k on M. 

Therefore, there exists a globally defined solution u(x,y) which has initial data 
u(x,y) = uq, d yi u(x,y) = u± on N, which vanishes to order k on M. Hence any 
other solution v(x, y) which takes on specified data on M up to fc-many derivatives 
may be changed by adding this solution u to it, without changing its initial data. 

Proof. We follow a construction that was used for the extension operators of section 
3. Let X3(£i v') be a Schwartz class function with support in the set j |r/ 1 2 < |£| 2 j C 
N. Its Fourier restriction to M, given by 



satisfies the constraint. While v may be nonzero on M, as may its derivatives, it 
is the case that for homogeneous polynomials pk+i(x", y") of degree k + 1, the 
function Pk+i {x" , y")v(x, y') on N vanishes on M to at least order k. Furthermore 
Pk+iv satisfies the constraint. Indeed, 



and differential operators do not affect the support. Set data Uo(x, y') = (pk+iv)(x, y') 
and ui — 0, and solve equation (1.2). Because this data satisfies the constraint, the 
solution u(x, y) is global. Because of the properties of the initial data, all x and 
y' derivatives of u(x, y) vanish on M. Because u\ = and u itself satisfies equa- 
tion (1.2), all y\ derivatives up to order k as well as any mixed derivatives also 
vanish. □ 





(Fpk+iv){(„rf) =Pk+i(-d£", TtV)X3 (£,»/). 
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5. A variant of a uniqueness theorem of Courant 

Courant (1962) gives a uniqueness result for the ultrahyperbolic equation with 
data posed on a hypersurface of mixed signature, which in our notation states that, 
among C 2 solutions, initial values of u(x,0,y') and d yi u(x, 0, y') prescribed in the 
set in the Cauchy hypersurface M given by 



^>-x°) 2 <« 2 , J>- % °) 2 < £ 2 



(5.1) 



1=1 1=2 

will determine a priori the values of the data on the larger set 



(x, y') G M : 



\ 



X>-%°) 2 < 



(5.2) 



1=2 



Furthermore the solution is determined uniquely in the space-time region 



(x,y) G 



l+d 2 



di 

t=i 



X>-y?) 2 <<* 



(5.3) 



Courant 's proof of this fact uses the Asgeirsson mean value theorem in a funda- 
mental way. 

The key implication from our point of view is that data on an arbitrarily small 
cylindrical subset of M, plus the stipulation of C 2 regularity, determine the data 
and indeed the solution on much larger sets of M and of space-time, respectively. 
In turn, knowledge of the data in a small cylinder determines the values of all of its 
derivatives on N — {(x, y) : y — 0} (if the data are smooth). This contrasts to the 
case discussed in section 4, in which it is shown that specification of a possibly large 
but finite number of derivatives does not lead to unique solutions, even when the 
constraint is imposed and the resulting solutions are globally defined and smooth. 

In this section we give a version of the above theorem of Courant, for data 
posed in ellipsoidal domains in the Cauchy hypersurface M, which are localized 
near the {y' — 0} coordinate axis (or any translate thereof). Our proof of this 
result is based on the Holmgren- John theorem (John 1982), and therefore remains 
true under perturbations to the equation. Thus it is a robust generalization of the 
Courant result, which being based on Asgeirsson's theorem is true only for precisely 
the ultrahyperbolic equation. 

Theorem 5.1. Let e > and define the ellipsoid Z E C M by 



{(x,y) : yi=0 , \x\' 



<1} 



< e < 1 . 



(5.4) 



A C 2 solution to (1.2) whose Cauchy data vanishes on Z E must necessarily vanish 
on the set 

D = {(x,y)eR dl+d2 : \x\ + \y\ < 1} 

and in particular its Cauchy data along with all derivatives must vanish on the 
subset of the Cauchy hypersurface given by {(x,y') £ M : \x\ + \y'\ < 1}. 
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Proof. Define R e (w) to be the cone over Z e with vertex v — (0, w\, w') G {(a;, y) : 
x = 0}. We will show that for any w — (wi,w') with |io| < 1 (namely the unit 
sphere in R d2 ), the region between the cone R s (w) and the ellipsoid Z e is a region 
of determinacy for the ultrahyperbolic equation. The closure of the envelope of such 
ellipsoidal cones includes the region D; in fact it is slightly larger. The result will 
follow accordingly. 

For a given R e (w), the Holmgren John theorem is based upon the construction 
of an analytic family of noncharacteristic hypersurfaces S\ with which to sweep 
the region between Z e and R e (w). Taking the case of the vertex v = (0,w) with 
w = ei := (1,0), define 

S x :={(x,y) : (1 - Vl f (\x\ 2 + Kf) = -A} 

with — 1 < A < 0. The normal to S x is N x = -2(x, (1 - Vi),y' /e 2 ) T , so that the 
characteristic form calculated on N x is 

Taking into account that (x, yi,y') <G S\ and solving for (1 — j/i) 2 , 

Recalling that A < (except in the limiting case S\ — » R e ) observe that this 
family of hyperboloids constitutes a noncharacteristic analytic family which sweeps 
the region between Z e and R s (e\). Thus the Holmgren- John uniqueness theorem 
applies, and this region is a region of determinacy for the ultrahyperbolic equation 
(1.2). 

We have already achieved the analogue of the statement (5.3) of Courant. 
Namely, given the values of a C 2 solution u(x, y) to (1.2) in the space-time ellipsoid 

W £ :={(x,y) : \x\ 2 + ^ < 1} , 

we may slice it with a hyperplane which contains the x-coordinate axes but which 
is otherwise arbitrarily oriented in y, to determine a possible Z e , which in turn 
determines the solution over the larger conical region R e with base Z e . All of these 
regions have been shown to be domains of determinacy. Their union contains the 
set D = {(x, y) : \x\ + \y\ < 1}. Therefore if a solution vanishes in W e it must also 
vanish in D. 

Returning to the problem of the domain of determinacy of the set Z e C M, we 
generalize the above construction to any w € R d2 with \w\ — 1. Let w — Re\ for 
ei = (1,0...), where R is an orthogonal matrix. Changing variables to z = Ry 
and using a symmetric matrix Q of signature (—, + ...), an analytic family of 
hyperboloids is given by 

S x (w):={(x,z) : \x\ 2 + ((z - e\),Q{z - e x )) = A} , 
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where the Euclidean inner product is given by (■,■)• The matrix Q is to be chosen 
so that the intersections of the hyperboloids S\(w) with the hypersurface M lie in 
Z £ , and sweeps it as A is varied. 

At this point we may assume without loss of generality that w = (u>i,w') = 
(toi, W2 1 . . . ), whereupon Q may be chosen such that 



for Q2 a 2 x 2 symmetric matrix with signature (—,+). Furthermore, the above 
rotation is then set to be 

_ (R 2 \ / cos(tf) sin(0) 

^0 I") ' 2 \-sm{6) cos(0) 

In y— coordinates the hypcrboloid family is expressed 

S x (w):={(x,y) : \x\ 2 + ((y - w), R T QR(y — w)) = A} , 

and the stipulation is that So (w) should intersect the hypersurface M in the original 
ellipsoid Z e . This imposes the condition that 

\x\ 2 + ((x,0,y'),R T QR(x,0,y')) := \x\ 2 + ((x,0,y'),B(x,0,y')) = \x\ 2 + ^\y'\ 2 , 

where B2 is the upper left-hand 2x2 block of the matrix B. Therefore one finds 
the matrix elements of B2 

e 2 -sin 2 (6») tan(0) 1 

611 ecos^e) ' 512 h2 -^' 

and furthermore, the 2x2 matrix Q 2 is 



_ / -1 tan(0)\ 
^ 2 ^tan(fl) ^ 



(5.5) 



where a = a(e,6) = (1 + (1 — e 2 )tan 2 (#)). Calculating the characteristic form on 
the hyperboloids S\(w), we compute the normal N\(w) as 

-±N x (w) = (x,Q(z- ei )) T . 

Noting that the characteristic form is invariant under rotations R as above, which 
leave the coordinate subspaces R^ 1 and R^ 2 invariant, we find that 

l N x (w) T (- Id ^ T ° )N x (w)=-\x\ 2 + ((z-e 1 ),Q 2 (z-e 1 )) • 



4 \ Id 2 xd 2 , 

This is evaluated on the hyperboloid S\(w), on which 

|x| 2 + ((z- ei ),Q(z- ei )) = A . 
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Solving for \x\ 2 , we find 

\N x ( W f (- 7 * xdl ^ J N x (w) = ((z ei ), [Q 2 + Q](z - ei )) - A . 
Specifically, the matrix [Q 2 + Q] is 

Using the form (5.5) for Q 2 , one calculates 

[o2 , 0] /tan 2 (9) £tan(0) \ 
W» + ^-^tan(tf) £+tan 2 (0)J " 

It is easily verified that this is positive definite. Recalling that A < in the def- 
inition of the analytic families of hyperboloids, it follows that S\(w) are all non- 
characteristic, and hence the Holmgren- John theorem applies, thus completing the 
argument. □ 
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